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SECTION  I 
INTRODUCTION 


Ever  since  Vander  Lugt  demonstrated  optical  matched 
filtering  in  1964  (Reference  1),  the  main  stream  of  optical 
pattern  recognition  has  been  correlation  type  image 
recognition.  This  field  of  research,  however,  has  had  limited 
success  in  practical  systems.  One  of  the  most  important 
reasons  was  the  lack  of  a  suitable  spatial  light  modulator 
(SLM).  The  development  of  the  Magneto-Optic  Device  (MOD)  has 
opened  up  new  possibilities.  This  report  addresses  the  issues 
of  the  feasibility  of  using  the  Magneto-Optic  Device  as  a 
spatial  light  modulator  for  correlation,  and  the  performance  of 
binary  filters  versus  linear  filters  in  terms  of  output  signal 
to  noise  ratio  (SNR)  and  storage  capacity.. 


R 


SECTION  II 
BINARY  CORRELATORS 


A  schematic  diagram  of  the  classic  Vander  Lugt  correlator 
is  shown  in  Figure  1.  The  real-time  implementation  of  this 
processor  requires  two  MODs .  The  first  is  used  to  record  the 
input  image  placed  at  plane  P^.  whereas  the  second  MOD. 
placed  at  plane  P  .  records  the  Fourier  transform  hologram  of 
the  reference  image.  The  correlation  output  is  obtained  at 
plane  P3> 

Since  the  MOD  is  a  binary  device.  the  effects  of 
thresholding  the  input  image  and  the  Fourier  transform  (FT)  of 
the  reference  image  have  to  be  evaluated.  The  performance 
measure  that  is  adopted  here  is  SNR.  the  signal  (correlation 
peak)  to  noise  (additive  noise  and  average  side-lobe)  ratio. 

The  first  investigation  involves  the  effect  of 

thresholding  the  input  image  on  the  SNR.  Let  f(i,j)  be  the 

image  to  be  recognized.  Assume  f(i,j)  to  be  a  discrete 

sequence  of  independent.  identically  distributed  Gaussian 

2 

random  variables  with  zero  mean  and  variance  .  The 
input  image  to  the  correlator  is  the  sum  of  f(i,j),  the  image 
to  be  recognized,  and  n(i,j),  the  additive  noise.  Here  it  is 
also  assumed  n(i.j)  to  be  an  independent,  discrete  sequence  of 

Gaussian  random  variables  with  zero  mean  and  variance 

2  ... 
on.  The  thresholded  image  is  defined  as 

1.  if  f(i,j)  +  n(i.j)  >0; 

f(i.j)  =  .  (l) 

-1,  if  otherwise. 

Assumed  for  now  is  a  joint  transform  architecture  where  both 
the  input  and  reference  images  are  recorded  on  MODs.  In  this 


2 


case,  the  reference  image  h(i.j)  is  obtained  by  thresholding 
f  ( i .  j ) : 


h(i.  j )  = 


1.  if  f (i.  j)  >  0; 

-1.  if  otherwise. 


(2) 


The  correlation  between  the  thresholded  input  and  reference 
images  is  then  given  by 


N  N 

g(i  '  »  3  1  )  =  Y  Y 

i  3 


f ( i. j )h( i  +  i  '  .  j  + j  '  ) 


The  signal  to  noise  ratio,  defined  as 


(3) 


SNR'  = 


E2  [g(0. 0)  ] 
var  [g(i ' . j '  ] . 


can  be  shown  to  be 


3’  =  0 


(4) 


SNR 


2 


N 


1 

2 


+ 


1  -1 
_  tan  ( SNR^n 

ir 


where 


SNR 


in 


(5) 


In  order  to  see  the  effect  of  thresholding  the  input 
image,  compare  SNR'  with  the  signal  to  noise  ratio  that  is 
obtained  if  the  input  and  reference  images 

thresholded.  It  can  be  shown  that  in  this  case 


are 


not 


SNR  = 


2  2 
N  SNR. 

m 


1  +  SNR 


The  two  SNRs  are  plotted  in  Figure  2  as  a  function  of 
2 

SNR^n>  As  can  be  readily  observed  from  the  plot,  the 
correlator  with  input  image  thresholded  is  only  marginally 
degraded  for  high  input  SNR. 

The  effect  of  thresholding  the  Fourier  transform  of  the 
reference  image  in  a  Vander  Lugt  correlation  in  which  the  FT 
hologram  is  recorded  on  an  MOD  is  explored.  (Since  the  Fourier 
transform  of  an  image  is  generally  complex,  what  is  meant  by 
thresholding  the  FT  is  that  the  real  part  of  the  FT  of  an  image 
pattern  is  hard-limited.)  Suppose  F_(u.v)  +■  jFT(u.v)  is 

the  FT  of  the  image  pattern  to  be  recognized.  The  computer 
generated  hologram  (CGH) ,  which  is  used  as  the  filter  for 
recognition,  is  then  formed  using  the  following  algorithm: 


H(u.v)  = 


+1.  if  Fr(u,v)  >  0; 
-1,  if  otherwise. 


The  correlation  output  g(x,y)  is  then  given  by 


g(x.y)  = 


V  u 

11 

-'-V.  u . 


||  Fr(u,v)  II  exp[  j2ir(ux  +  vy)]  dudv 


1  “1 


-v,  u 


f  1  f  1 

/  /  F.(u.v)  H(u,v)  exp[]2ir(ux  +  vy)  ]  dudv 

v ,  u. 


1  1 


The  first  term  on  the  right-hand  side  of  Equation  8  is  the 
Fourier  transform  of  a  real,  positive  quantity  and  will. 


Rtf 


therefore,  give  a  strong  correlation  peak  at  (x.y)  =  (0,0). 
The  expected  value  of  the  second  term  is  zero.  Hence,  a  binary 
CGH  of  this  type  can  be  used  for  recognition.  Note  that  the 
usage  of  such  an  on-axis  CGH  is  possible  only  because  the 
bipolar  nature  of  the  CGH  eliminates  the  necessity  of  bias 
removal . 


The  loss  in  SNR  that  results  from  the  use  of  this 

suboptimum  filter  is  estimated.  Let  the  image  f(x.y)  to  be 

recognized  be  a  Gaussian  white  process  with  zero  mean  and 
2 

variance  .  Suppose  f(x,y)  is  contaminated  by  an 

additive  Gaussian  white  noise  n(x.y).  independent  of  f(x,y), 

2 

with  zero  mean  and  variance  a  .  For  a  »  ac,  it 

n  n  f 

can  be  shown  (Reference  2)  that: 


SNR  = 


where  SBP  is  the  space-bandwidth  product.  Recall  that  the  SNR 
of  the  conventional  matched  filter  is  the  product  of  the  SBP 
and  the  input  SNR.  Thus,  the  SNR  is  only  reduced  by  a  factor 
of  ir  when  a  quantized  hologram  replaces  a  matched  filter.  It 
is  difficult  to  find  a  closed  form  expression  for  the  SNR  for 
sufficiently  low  noise  (o£  >>  a  ) .  However.  numerical 
evaluations  for  low  SBP  indicate  that  the  performance  of  the 
quantized  FT  filter,  while  always  slightly  inferior  to  that  of 
the  matched  filter,  approaches  the  performance  of  the  matched 
filter  as  the  input  SNR  becomes  large  (Figure  3).  Experimental 
evidence  also  supports  this  result  for  higher  SBP.  The  results 
of  two  optical  correlation  experiments  using  the  MOD  programmed 
with  quantized  FT  holograms  are  shown  in  Figures  4  and  5. 


The  MOD  can  also  be  used  in  the  Fourier  plane  of  a 
two-dimensional  holographic  acousto-optic  image  correlator. 
Figure  6  shows  the  schematic  diagram  of  the  corresponding 
architecture.  The  cross-correlation  is  obtained  by  taking  the 
inverse  Fourier  transform  of  the  product  of  the  Fourier 
transforms  of  the  input  image  and  the  reference  in  the  abscissa 
while  performing  a  shift  and  integrate  type  correlation  in  the 
ordinate.  The  procedure  for  producing  the  correlation  can  be 
described  mathematically  by  the  following  equation: 

N 

G(u.y)  =  £  F(u-  nd  +  7)  H(u.nd)  (10a) 

n 

and 


g(x.y) 


/a<  “• 


y)ej2irUX  dx, 


(10b) 


where  F(u,nd  +  d)  and  H(u,nd)  are  the  one-dimensional  Fourier 
transforms  of  the  input  image  and  the  reference,  respectively. 
As  it  is  apparent  from  the  equation,  the  algorithm  for 
generating  the  hologram  is  different  from  the  one  described  by 
Equation  7.  Since  the  correlation  operation  of  this 
architecture  involves  space  integration  (implemented  by  the 
one-dimensional  Fourier  transform)  in  one  dimension  and  time 
integration  (implemented  by  the  time  delay  integrating  mode  of 
the  charge  coupled  device  (CCD))  in  the  other,  the  hologram  is 
generated  by  quantizing  the  real  part  of  only  the  one¬ 
dimensional  Fourier  transform  of  the  reference  image;  that  is: 


H(u , nd ) 


+  1. 

-1. 


if  F  (u  ,  nd  )  >_  0; 
if  F  (u  .  nd  )  <  0 
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SECTION  III 

BINARY  ROTATION  INVARIANT  FILTERS 


Mathematically.  any  pattern  f(r,0)  can  be  decomposed 
into  an  infinite  sum  of  its  circular  harmonics. 


f(r.e)  =  X  f_(r)e  jne 


(12a) 


where 


fn(D  = 


J f(r- 


o)e5n6de>. 


(12b) 


Thus,  if  the  pattern  is  correlated  with  any  component  of  its 

-  i  n0 

circular  harmonics  fn(r)e  .  the  magnitude  of  the 
output  is  rotation  invariant  (Reference  3).  Furthermore,  if 
any  component  of  the  circular  harmonics  is  recorded  on  a 
Fourier  transform  hologram,  the  resulting  filter  is  also  shift 
invariant . 

Bipolar  filters  can  be  formed  by  using  an  algorithm 
similar  to  that  described  by  Equation  7,  namely: 


H(s.w)  = 


+1.  if  Re [ F  (s ,w) ]  >  0; 
n 

-1,  otherwise. 


(13) 


where  F  (s.w)  is  the  Fourier  transform  of  the  nth  order 
n 

circular  harmonics  of  the  input  image.  Computer  generated 
bipolar  filters  were  obtained  from  the  first  two  orders  of 
circular  harmonics  of  the  letter  A.  Numerical  correlations  of 
the  letters  A.  B.  C,  and  D  and  the  bipolar  CGHs  were 
evaluated.  The  result  are  presented  in  Figure  8.  The 
intra-class  recognition  using  the  zeroth  order  circular 


■  fc’***  •  »  k  *  j,  *  , 


harmonics  is  remarkable.  However,  the  inter-class 
discrimination  performance  is  not  satisf actory .  In  order  to 
obtain  good  discrimination  performance,  cross-correlations 
using  binary  filters  generated  from  different  circular 
harmonics  have  to  be  evaluated.  A  compromise  between  the  inter- 
and  intra-class  discrimination  performance  will  determine  which 
binary  filter  is  to  be  used.  The  result  of  the  optical 
implementation  of  the  cross-correlation  of  the  letter  A  and  the 
binary  Fourier  transform  filter,  recorded  on  a  128  by  128  MOD, 
of  the  A 1  s  zeroth  order  circular  harmonics  is  presented  in 
Figure  9.  Figure  9(a)  is  the  image  of  the  rotation-invariant 
binary  filter,  whereas  Figure  9(b)  is  the  reconstruction  of  the 
binary  filter.  Note  that  both  of  them  are  circularly 
symmetric.  Figure  9(c)  is  the  correlation  output.  Note  that 
the  peak  is  very  sharp  and  is  located  at  the  center  of  the  dc 
pass  of  the  A,  just  as  the  simulation  predicted.  Rotation  of 
the  input  object  did  not  cause  the  peak  to  diminish  in 
amplitude  as  expected. 


SECTION  IV 

GENERALIZED  BINARY  FILTERS 


This  section  examines  the  possibility  of  designing  linear 
bipolar  filters  for  recognizing  multiple  objects.  One 
immediate  extension  of  such  a  design,  if  it  can  be  done,  is  to 
incorporate  information  corresponding  to  multiple  versions  of  a 
certain  object  in  a  filter  to  achieve  invariance  (e.g., 
rotation  and/or  scale).  Two  schemes  that  may  accomplish  this 
goal  are  investigated  as  described  below. 

The  first  scheme  is  a  derivative  of  the  well  known 
2-category  perceptron  invented  in  the  1960’s  (Reference  4). 
Basically.  a  2-category  linear  perceptron  is  a  machine 
consisting  of  a  weight  pattern  w  trained  by  a  certain  algorithm 
utilizing  the  set  of  all  input  patterns  {f^}  such  that 
the  output 


N2 

g  =  sgn 

.  i.  j  J 


+1  if  f  belongs  to  Y; 
m 

-1  otherwise. 


(14) 


Attention  is  here  confined  to  one  of  the  simplest  perceptron 

training  algorithms,  which  is  as  follows:  the  set  of  input 

patterns  is  arranged  as  a  repetitive  sequence  of  patterns, 
12  Ml 

i.e.,  {f_  .  £ .  f  .  f.  ...}.  The  weight  vector 

is  trained  by  the  elements  from  the  sequence  one  at  each  time. 


K  -+■  1  Jc  jc  k  1 

w  =  w  +af;w  =  some  initial  vector. 


(15a) 


i'iw 


where 


+  1  if 

k 

k 

k  k 

g 

=  sgn  [w 

1  ]  =  -1  and  f_  e  V; 

-1  if 

g* 

=  +1  and  f k 

$ 

(15b) 

0  otherwise. 


A  solution  vector  is  found  when  the  right  response  is  obtained 
for  each  input  vector. 

A  bipolar  perceptron  can  be  defined  by  modifying  the 
linear  perceptron  simply  by  replacing  w  by  a  bipolar  weight 
pattern.  The  training  algorithm  can  also  be  modified  to 


„k  +  l  _  s g n  [w^  +  a^f*]  (16) 

The  result  of  the  rate  of  convergence  of  a  binary 

perceptron  trained  by  using  the  above  scheme  is  shown  in  Figure 

4  . 

10.  An  ordinate  reading  of  10  iterations  means  the  binary 

4 

perceptron  did  not  converge  within  10  steps.  Even  though  a 
convergence  proof  does  not  exist,  preliminary  computer 
simulations  indicate  that  a  binary  solution  weight  pattern  can 
be  obtained  provided  the  dimensionality  N  is  high  enough.  The 
statistical  capacity  with  N  =  16  is  shown  in  Figure  11.  Note 
that  the  capacity  curve  starts  to  roll  off  at  approximately  M. 
the  number  of  stored  vectors  equal  to  2N/3,  whereas  the 
capacity  curve  of  a  linear  perceptron  (not  shown)  does  not  roll 
off  until  M  =  N.  This  can  be  explained  by  the  fact  that  the 
solution  region  of  a  linear  perceptron  may  not  contain  a  binary 
point.  In  such  a  case,  a  binary  vector  does  not  exist  which 
classifies  all  input  vectors  correctly. 
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The  second  scheme  that  can  be  used  to  recognize  multiple 
patterns  is  described  below.  For  ease  of  comparison,  the 
ordinary  (non-bipolar )  linear  filter  of  the  second  scheme, 
designed  for  the  desired  purpose,  is  first  analyzed.  Consider 
the  following  algorithm:  the  operation  to  be  performed  is 

given  as 


g  =  sgn  £  Mi. 5) 


+1  if  f  belongs  to 
-1  otherwise. 


(17) 


where  ¥  is  the  class  of  objects  to  be  recognized.  The  filter 
h(i,j)  is  generated  by  forming  a  weighted  sum  of  all  the  input 
patterns,  i.e.. 


ha.j)  -  £  am,  £„,  (i.j) 


(18a) 


where 


+  1.  if  f  ,  belongs  to 
m 

-1,  otherwise. 


(18b) 


To  see  that  the  above  scheme  is  capable  of  recognizing  multiple 
objects,  rewrite  Equation  17  as 


yc2m  a.)) .  2  2  <i9> 

i.j  m '  ><m 


The  right-hand  side  of  Equation  19  is  composed  of  two  terms, 
namely,  the  signal  (first)  term  and  the  noise  (second)  term. 
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Provided  that  the  signal  term  is  sufficiently  large  compared  to 
the  noise  term,  the  correct  response  is  expected. 

The  next  task  is  to  obtain  a  theoretical  estimate  of  the 

capacity  of  the  filter.  The  capacity  is  defined  to  be  the 

number  of  vectors  M  that  can  be  stored  in  a  filter  with 

c 

vanishingly  small  probability  of  error  for  sufficiently  large 

dimensionality  N.  Assuming  f  (i.j)  to  be  a  discrete  sequence 

m 

of  bipolar-valued  independent  random  variables,  i.e. 


P[f  (i.j)  =  1]  =  P[f  (i.j)  =  1]  =  1/2; 
m  m 

P[fm(i'j)  1  fm-(i''j,)]  = 


(20) 


it  can  be  shown  that 


Mc  = 


8  in  N 


(21) 


Shift  invariance  can  be  incorporated  into  this  system  by 
modifying  the  operation  to  be 


g(i'.j')  =  sgn  £  h( i . j ) f  ( i  +  i  '  .  j  +  j') 


(22) 


+1  if  f  belongs  to  V; 
m 

-1  otherwise. 


In  such  case,  the  capacity  Mc  can  be  shown  to  be 


Mc  = 


16  an  N 


(23) 


The  digitally  computed  correlations  of  random  sequences 
(statistics  given  by  Equation  20)  and  the  linear  weighted-sum 
filter  with  different  number  of  vectors  stored  are  shown  in 
Figure  12.  It  can  be  seen  from  the  simulated  results  that  as 
the  number  of  vectors  stored  in  a  filter  increases,  so  do  the 
magnitudes  of  the  sidelobes.  This  phenomenon  accounts  for  the 
limit  of  the  number  of  vectors  that  can  be  stored  in  a 
weighted-sum  filter.  The  histogram  of  capacity  with 
dimensionality  N  equal  128  is  shown  in  Figure  13.  The  mean  of 
the  histogram  agrees  with  the  theoretical  result. 

For  the  bipolar  analog  of  the  above  scheme,  consider  the 
same  operation  given  by  Equation  17.  In  this  case,  the  filter 
is  generated  by  bipolarizing  the  filter  given  by  Equation  18. 
It  is  given  as  follows. 


(24) 


h(i.  j)  -  sgn  £  am,  f m,(i.  j) 

m 1 


Using  the  same  statistics  for  the  stored  vectors,  it  can  be 
shown  that  the  capacity  of  a  bipolar  weighted  sum  filter  is 
only  reduced  by  a  factor  of  ir/2.  Thus 


Mb  =  — 


w  16  In  N 


for  filters  with  shift  invariance  incorporated,  and 


M  =  2  N 

b  if  8  In  N 


(25) 


(26) 


for  filters  without.  Computer  simulated  correlations  of  random 
vectors  and  the  bipolar  filter  formed  by  the  above  algorithm 
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are  shown  in  Figure  14.  Note  that  the  magnitude  of  the  peak 
decays  as  the  number  of  vectors  stored  increases  in  this  case. 
This  phenomenon  is  understandable,  since  as  M  increases,  the 
relative  information  of  each  vector  stored  in  the  filter 
decreases.  The  histogram  of  capacity  of  the  corresponding 
bipolar  filter  is  shown  in  Figure  15.  The  mean  of  the 
histogram  agrees  with  the  theoretical  result. 


SECTION  V 

MULTIPLE  BITS  CORRELATORS 


Conventional  matched  filtering  can  be  implemented  by 
utilizing  a  binary  device  such  as  the  MOD.  Let  the  input  image 
be  f(i,j).  The  reference  image  h(i,j)  can  be  approximated  by  a 
weighted  sum  of  binary  numbers  as  the  following 


h( i . j )  »  £  hk  (i, j)2*  .  (27) 

k=0 

The  correlation  output  is  thus  given  by 


The  term  included  in  the  brackets  is  simply  the 

correlation  of  the  input  image  and  a  binary  filter.  Thus,  the 
final  output  will  be  the  sum  of  all  the  K  weighted 

correlations,  each  being  computed  separately. 

One  way  to  implement  this  scheme  is  to  have  a  setup  of 
multiple  MODs,  each  performing  a  bit-slice  correlation  with  the 
input  image.  All  the  correlations  are  weighted  properly  and 
are  then  summed  up  interf erometrically .  This  implementation 
has  the  drawback  of  having  to  use  multiple  MODs.  If  the  Vander 
Lugt  correlation  is  to  be  computed,  a  total  of  2K  MODs  is 
required  (K  MODs  for  the  real  part  of  the  FT  of  h(i,j)  and 
another  K  MOD  for  the  imaginary  part). 
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Another  way  to  implement  the  bit-slice  scheme  is  to 
compute  the  bit-slice  correlation  in  time  by  updating  the  MODs 
(two  are  required  in  this  case)  with  different  bit  planes  of 
the  reference  image.  The  final  output  is  obtained  by 
accumulating  all  the  properly  weighted  correlations  at  the 
output  detector.  The  tradeoff  here  is  that  it  requires  K  times 
as  long  to  perform  a  conventional  correlation  as  it  takes  to 
compute  a  bit-slice  correlation.  The  result  of  the  computer 
simulation  of  the  second  method  of  implementation  is  shown  in 
Figure  16.  Note  that  the  SNR  of  the  final  output  differs  from 
the  first  bit-slice  correlation  by  a  factor  of  approximately 
three.  In  view  of  this  and  similar  results,  and  in  addition, 
the  knowledge  gained  from  the  analysis  of  the  previous 
sections,  it  is  believed  that  the  gain  in  the  output 
performance  by  computing  bit-slice  correlation  does  not 
rightfully  compensate  the  increase  in  cost,  difficulty  and 
complexity. 
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SECTION  VI 
CONCLUSION 
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This  report  has  shown  that  the  performance,  in  terms  of 
output  SNR  and  storage  capacity,  of  binary  correlation  filters 
does  not  severely  degrade  when  compared  to  conventional  correla¬ 
tion  filters.  As  a  matter  of  fact,  bipolar  spatial  light 
modulators  (SLMs)  are  free  of  the  problems  of  limited  dynamic 
range  and  non-1 inear i ty.  However,  binary  filters  are  not 
without  other  problems  besides  those  mentioned  earlier.  Many 
issues  still  remain  to  be  researched,  for  instance,  how  a  binary 
filter  is  optimized.  How  a  binary  perceptron  is  trained  to 
recognize  multiple  objects  is  another  example.  For  a  certain 
set  of  data,  a  binary  filter  which  will  make  the  right  decision 
in  recognition  may  not  exist  at  all.  It  is,  nevertheless,  the 
authors'  expectation  that  the  probability  of  having  such  a 
binary  filter  approaches  one  for  sufficiently  high 
dimensionality. 
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(a)  INPUT 


(b)  HOLOGRAM 


(c)  CORRELATION  OUTPUT 


7.  Experimental  Demonstration  of  the  Holographic 
Acoustic  Correlator  with  a  Programmable 
Magneto-Optic  Device 
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CONVERGENCE  RATES  VS.  NUMBER  OF  VECTORS 


TOTAL  NUMBER 
OF  ITERATIONS 
FOR  PECEPTRON 
TO  CONVERGE 


NUMBER  OF  VECTORS  (  dim  =  63) 


Figure  10.  The  Rate  of  Convergence  of  a  Bipolar  Perceptron 
(The  Rate  of  Convergence  of  an  Ordinary  Perceptron  of 
the  Same  Data  Set  is  Also  Plotted  for  Comparison.  The 
Circles  on  the  Plot  are  the  Bipolar  Vectors 
Generated  by  Taking  the  sgn  []  of  the 
Ordinary  Perceptron  that  Dichotomize 
the  Data  Set  Correctly) 
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Icl  CORRELATION  OUTPUT  (N  =  128,  M  =  4) 
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Id)  CORRELATION  OUTPUT  IN  -  128,  M  -  6) 


Figure  12.  The  Simulated  Correlations  of  the  Linear 
Weighted-Sura  Filters  of  Different  Number  of  Stored 
Vectors  and  Some  of  Their  Stored  Vectors  (The 
Threshold  Level  for  Recognition  is  One 
Half  of  the  Energy  of  Each  Vector) 
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Figure  13.  The  Histogram  of  Storage  Capacity  Based 

100  Trials 
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NUMBER  OF  PIXELS 


(cl  CORRELATION  WITH  BIPOLAR  FILTER  (M  -  31 
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Id)  CORRELATION  WITH  BIPOLAR  FILTER  (M  -  41 


Figure  14.  Computer  Simulated  Correlations  of  Bipolar 
Weighted-Sum  Filters  of  Different  Number  of  Stored 
Vectors  and  Some  of  Their  Stored  Vectors 
(The  Threshold  Level  for  Recognition  is 
One-Half  of  the  Energy  of  Each  Vector) 


Numerical  Evaluation  of  Bit-Slice  Correlations 
at  Different  Time  Integration  Intervals 
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